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Abstract 



In this paper we consider the Cahn-Hilhard equation endowed with Wentzell 
boundary condition which is a model of phase separation in a binary mixture con- 
tained in a bounded domain with permeable wall. Under the assumption that the 
nonlinearity is analytic with respect to unknown dependent function, we prove the 
convergence of a global solution to an equilibrium as time goes to infinity by means 
of a suitable Lojasiewicz-Simon type inequality with boundary term. Estimates of 
convergence rate are also provided. 

Keywords: Cahn-Hilliard equation, Wentzell boundary conditions, Lojasiewicz- 
Simon inequality, convergence to equilibrium. 

1 Introduction 

This paper is concerned with the asymptotic behavior of the global solution to the fol- 
lowing Cahn-Hilliard equation 



in [0, T] X Q 



(1.1) 



II = -Au + f{u), 



in [0, T] X Q, 



(1.2) 
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subject to Wentzell boundary condition 

Afi + bd^fi + c/i = 0, on [0, T] x T, (1.3) 
the variational boundary condition 

-aA\iu + d^u + (3u= ^, on [0,T]xr, (1.4) 

and initial datum 

u{0,x)=iIjo, in fl. (1.5) 

In above, < T < oo, is a bounded domain in (n = 2, 3) with smooth boundary T. 
a,P,b,c are positive constants. Ay is the Laplace-Beltrami operator on F, and z/ is the 
outward normal direction to the boundary. 

The Cahn-Hilliard equation arises from the study of spinodal decomposition of binary 
mixtures that appears, for example, in cooling process of alloys, glass or polymer mixtures 
(see [1,12,20,28] and the references cited therein), /i is called chemical potential in the lit- 
erature. The classical Cahn-Hilliard equation is equipped with the following homogeneous 
Neumann boundary conditions 

dufi = 0, t>o, xe r, (1.6) 
d^u = 0, t>o, xe r. (1.7) 

Boundary (II. 6p has a clear physical meaning: there cannot be any exchange of the mixture 
constituents through the boundary T which implies that the total mass udx is conversed 
for all time. The boundary condition (ll.7p is usually called variational boundary condition 
which together with 01.61) result in decreasing of the following bulk free energy 

Eb{u) = j Q|Vm|2 + F(m)^ dx, (1.8) 

where F{s) = f{z)dz. A typical example in physics for potential F is the so-called 
'double-well' potential F{u) = |(m^ - 1)^. 

For the equations (11.11) (II. 2p subject to boundary conditions HIM (11.71) and initial 
datum (11.51) . extensive study has been made. We refer e.g., to [3,7,19,25,28,36] and the 
references cited therein. In particular, convergence to equilibrium for the global solution 
in higher space dimension case was proved in [25]. 

Recently, a new model has been derived when the effective interaction between the 
wall (i.e., the boundary F) and two mixture components are short-ranged (see Kenzler et 
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al. [12]). In such a situation, it is pointed out in [12] that, the following surface energy 
functional 



Es{u) = J [f r ll^r + f « - hsu) dS, (1.9) 

with V|| being the covariant gradient operator on F (see e.g. [17]), should be added to the 
bulk free energy Ef){u) to form a total free energy functional 

E{u)=Eb{u) + E,{u). (1.10) 

In above, o-g > 0, gs > 0, hg are given constants. Together with the no-flux boundary 
fll.Gp condition, the following dynamical boundary condition is posed in order that the 
total energy E{u) is decreasing with respect to time: 

(TsA||M — d^u — QgU + hg = —Ut, t > 0, X eT. (1-11) 

J- s 

We refer to [2,17,23,24,29] for extensive study for system (ll.ipp.2p with boundary 
conditions (11.60(11.111) and initial datum 01.51) . In particular, Wu & Zheng [29] proved 
the convergence to equilibrium for a global solution as time goes to infinity by deriving a 
new type of Lojasiewicz-Simon inequality with boundary term (see also [2] for a different 
proof). 

Based on the above model, in a quite recent article by Gal [4], the author proposed 
( ll.ip - p.5p as a variation model which describes phase separation in a binary mixture 
confined to a bounded region fl with porous walls. Instead of the no-fiux boundary 
condition (II. 6p . the Wentzell boundary condition (II. 3p is derived from mass conservation 
laws that include an external mass source (energy density) on boundary T. This may be 
realized, for example, by an appropriate choice of the surface material of the wall, i.e., 
the wall r may be replaced by a penetrable permeable membrane (ref. [4]). Then, the 
variational boundary condition ( ll.4p is introduced in order that the system ( ll.ll )- (11.5p 
tends to minimize its total energy 



Namely, 



E{u)= I [-\\/u\' + F{u)]dx+ I [-\Viiu\ +^u' ] dS. (1.12) 



^E{u{t)) = - [ \Vfi\^dx -if fi^dS < 0. (1.13) 
dt Jq J-p 

For more intensive discussions, we refer to [4-6]. 

In [4], the existence and uniqueness of global solution to problem (|l.ip - (ll.5p has been 
proved by adapting the approach in [24]. Later, in [5], the same author studied the 
problem in a further way that he obtained the existence and uniqueness of a global 



3 



solution to the problem under more general assumptions than those in [4]. He showed 
that the global solution defines a semiflow on certain function spaces and also proved the 
existence of an exponential attractor with finite dimension. 

Then a natural question is: whether the global solution of system (ll.ip - (ll.5p will 
converge to an equilibrium as time goes to infinity? This is just the main goal of this 
paper. Moreover, we shall provide estimates for the rate of the convergence (in higher 
order norm). 

Remark 1.1. Without loss of generality, in the following text, we set positive constants 
b, c, a, P to be 1. In this paper, we simply use \\ ■ \\ for the norm on L'^{Q) and equip H^{Q) 
with the equivalent norm 

\\u\\m{Q)=(^J \Vu\'^dx + J u^dS^ . (1.14) 

Before stating our main result, first we make some assumptions on nonlinearity /. 
(Fl) f{s) is analytic in s G M. 
(F2) 

\fis)\<Cil + \s\n, VsGR, 

where C > 0, p > and p G (0, 5) for n = 3. 
(F3) 

liminf f'{s) > 0. 

>oo 

Remark 1.2. Assumption (Fl) is made so that we are able to derive an extended 
Lojasiewicz-Simon inequality to prove our convergence result. Assumption (F2) implies 
that the nonlinear term has a subcritical growth. Assumption (F3) is some kind of dis- 
sipative condition. (F3) is supposed in [4, 5] to obtain the existence and uniqueness of 
global solution to the evolution problem (ll.ip - (ll.5l) . Moreover, (F3) together with (F2) 
enable us to prove the existence result for stationary problem (ll.lSp by variational method 
(see Section 3). It's easy to check that the nonlinearity f{u) = u^ — u corresponding to the 



most important physical potential F{u) = |('U^ — 1)^ satisfies all the assumptions stated 
above. 



Let V be the Hilbert space which, as introduced in [4], is the completion of C^{Q) 
with the following inner product and the associated norm: 



{u,v)v 



/ Vu-Vvdx+ / {V\\u-V\\v + uv)dS, yu,veV. (1.15) 
Jn Jr 



The main result of this paper is as follows. 
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Theorem 1.1. Let (F1)-(F3) he satisfied. For any initial datum Uq G V, the solution 
u{t,x) to problem ^1.1\) - (T7^) converges to a certain equilibrium ip{x) in the topology of 
H^{VL) n H^{T) as time goes to infinity, i.e., 



lim {\\u{t, ■) - i)\\Hi(^n) + \\u{t, ■) - V'||//3(r)) = 0. 



(1.16) 



Moreover, we have the following estimate for the rate of convergence: 



u - nw^in) + \\u - nn^T) + \\ut\\v < C{1 + t)-'l^^-^'\ 



t > 1. 



(1.17) 



Here, C > 0, il){x) is an equilibrium to problem ( fj.jj j-( f773]i . i.e., a solution to the following 
nonlinear boundary value problem: 



Before giving the detailed proof of Theorem 1.1, let's first recall some related results 
in the literature. The study of asymptotic behavior of solutions to nonlinear dissipative 
evolution equations has attracted a lot of interests of many mathematicians for a long 
period of time. Unlike in 1-d case (see [16,35]), the situation in higher space dimension 
case can be quite complicated. On one hand, the topology of the set of stationary solutions 
can be non-trivial and may form a continuum. On the other hand, a counterexample has 
been given in [22] for a semilinear parabolic equation saying that even the nonlinear 
term being C°° cannot ensure the convergence to a single equilibrium (see also [21]). In 
1983 Simon in [26] proved that for a semilinear parabolic equation if the nonlinearity 
is analytic in unknown function u, then convergence to equilibrium for bounded global 
solutions holds. His idea relies on generalization of the Lojasiewicz inequality (see [13-15]) 
for analytic functions defined in finite dimensional space M™. Since then, Simon's idea has 
been applied to prove convergence results for many evolution equations, see e.g., [8-11,25] 
and the references cited therein. To the best of our knowledge, most previous work 
are concerned with evolution equations subject to homogeneous Dirichlet or Neumann 
boundary conditions. 

Our problem p.ll) - (|1.5p has the following features. The first boundary condition (11.31) 
is Wentzell boundary condition which involves the time derivative of u; the second bound- 
ary condition (11. 4p for u has a mixed type since it also involves the chemical potential 
11. It turns out that for the corresponding elliptic operator, it yields a non-homogeneous 




Az/^ + f{i)) = 0, X G 

A||^ + 9^V^ + V^ = 0, xer. 



(1.18) 



and 9 G (0, \) is a constant depending on il){x). 
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boundary condition. The Lojasiewicz-Simon inequality for homogeneous boundary con- 
ditions in the hterature fails to apply. As a result, a non-trivial modification is required 
to treat the present problem. We succeed in deriving an extended Lojasiewicz-Simon 
type inequality involving boundary term, with which we able to show the convergence 
result (for other applications, see [29-31,33]). Besides, by delicate energy estimates and 
constructing proper differential inequalities, we are able to obtain the estimates for the 
convergence rate (in higher order norm). This in some sense improves the previous re- 
sult in the literature (see for instance [8, Theorem 1.1]) and can apply to other evolution 
equations (ref. [31,32]). 

The rest part of this paper is organized as follows: In Section 2 we introduce the 
functional settings and present some known results on existence and uniqueness of global 
solution and uniform compactness obtained in [4,5]. In Section3 we study the stationary 
problem. Section 4 is devoted to prove an extended Lojasiewicz-Simon inequality with 
boundary term. In the final Section 5 we give the detailed proof of Theorem 1.1. 



2 Preliminaries 

We shall use the functional settings introduced in [4,5]. 

For u G C{Q), we identify u with the vector U = {u\^,u\^) G C{Q) x C(r). We define 

H = L'^{Q) © i^^(r) to be the completion of C{Q) with respect to the following norm 

1 

h\\H= (hf + Iklli2(r))'. (2.1) 

For any g eH, consider the elliptic boundary value problem 

—Am = fif, in fi, , , 

(2.2) 

d^u + u = g, on F. 
We can associate it with the following bilinear form on H^{Q): 



a{u, v) = Vm ■ Vudx + / uvdS, (2.3) 
Jq Jv 

for all M, ?; G H^iVt). Then it defines a strictly positive self-adjoint unbounded operator 
A ■ D{A) = {u G H\n)\Au eH}^n such that 

{Au,v)n = a{u,v), y u e D{A), v e H\n). (2.4) 

Then by Lax-Milgram theorem, it follows that the operator A is a bijection from D{A) into 
H and A~^ : H ^ H is a linear, self-adjoint and compact operator on H (see [5, Section 
2] or [4, Section 4]). In other words, for any g eH, A~^g is the unique solution to (12.21) . 
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We can then consider the weak energy space X endowed with the following norm, 

\\v\\l = \\A-^/\\\^^ = {A~\, v)n, yven. (2.5) 

It follows that 

{u,v)^ = {u,A-\)^, yueH\n),veX. (2.6) 
In particular, for all f G X and u = A~^v 

||f llx = {A^'^v,v)'j-c = {u,Au)n = a{u,u). (2.7) 

For more detailed discussions, we refer to [4,5]. 

The existence and uniqueness of global solution to p.ip -( fT75l l has been obtained in 
[4,5]. The results in [5, Section 3,4] and [4, Section 4] in particular imply that 

Theorem 2.1. Let (F1)-(F3) he satisfied. For any initial datum uq G V, problem 
f ll.ip - p.5l) admits a unique global solution u{t,x) which defines a global semiflow on V. 
Moreover, u{t,x) belongs to for t > 0. 

The total free energy 

E{u) = j Q|Vm|2 + F(m)^ + ^ Q|V||m|' + ^m2^ rfS. (2.8) 

where F{s) = f{z)dz, serves as a Lyapunov functional for problem (ll.ip - fll.5p . In 
other words, for the smooth solution u to problem (ll.ip - (ll.5l) . we have 

d 



,E{u)+ / |V/ir+ / \i2\'dS = 0. (2.9) 
dt Jq Jr 

Uniform bounds for the solution which yield the relative compactness in H^{Q)nH^(T) 
can be seen from [5, Proposition 3.3, Theorem 3.5], here we state the result without proof. 

Lemma 2.1. Let (F1)-(F3) hold and 7 G [0, 1/2). Then, for any initial datum Uq G V, 

the solution of fll.ip - (ll.5p satisfies the following dissipative estimates, namely, for any 
6 > 0, there hold 

Mt)\\ls+,^n) + lkWII|3+.(r) <Cs, t>6>0, (2.10) 
and (ref [5, (3.34)]), 

WMmm + \\ut\\m{v) <C5, t>5>0, (2.11) 
where Cs > depends only on \\uq\\v and S. 
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For any initial datum Mq ^ V, the tu-limit set of Uq is defined as follows: 



uj{uo) = {i'ix) I 3 {tn} such that u(tn,x) — > ipi^x) G V, as t„, — +00}. 
Then we have 

Lemma 2.2. For an?/ uq G l^, t/te u-limit set of uq is a compact connected subset in 
H^{n) n H^iT). Furthermore, 

(i) uj{uo) is invariant under the nonlinear semigroup S(t) defined by the solution u{x,t), 
i.e, S{t)u!{uo) = uj{uo) for all t > 0. 

(a) E{u) is constant onuj{uQ). Moreover, uj{uq) consists of equilibria. 

Proof. Since our system has a continuous Lyapunov functional E{u), the conclusion of 
the present lemma follows from Lemma 12.11 and the well-known results in the dynamical 
system (e.g. [28, Lemma LLl]). Thus, the lemma is proved. □ 



3 Stationary Problem 

In this section we study the stationary problem. The stationary problem corresponding 
to (iri|)-([L5]) is 

f 

Afi = 0, X eQ, 
-A^/j + f{^) = J1, xen, 
duji + /I = 0, a; G r, 
-Aii'ip + d^ip + ip = ll, X G r. 

Then it immediately follows that Jl = and the stationary problem is reduced to p.l8p . 

Lemma 3.1. Let (F1)-(F3) be satisfied. Suppose that G H^{Q) n H^{T) satisfies 
^1.18\) . Then ip is a critical point of the functional E{u) over V . Conversely, if ip E V is 
a critical point of E{u), then ip G C°° and it is a classical solution to problem ^1.18\) . 

Proof. The proof is similar to [29, Lemma 2.1]. The C°° regularity for solution ip follows 
from the elliptic regularity for (11.181) (see e.g., [17, Corollary A.l, Lemma A. 2]) and a 
bootstrap argument. □ 

Lemma 3.2. Let (F1)-(F3) be satisfied. The functional E{u) has at least a minimizer 
V E V such that 

E{v) = inf E{u). (3.2) 
In other words, problem ^.18\} admits at least a classical solution. 
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Proof. From assumption (F3), there exists 5o > such that 



liminf /'(s) > (5o. (3.3) 



Then there exists A^i = Ni{5o) > such that 

1 
2 



f'is) > l6o, \s\ > N,. (3.4) 



It follows that 

liminf /(s) > 1, liminf /(s) < -1. (3.5) 

s— >+oo s— >— oo 

Since F'{s) = /(s), then we can deduce from (13.51) that 

liminf F(s) > 1. (3.6) 

s|— »+oo 

Therefore, there exists iV2 > such that 

Fi-s) > 0, |s| > N2. (3.7) 

This indicates that 

F{u)dx= / F{u)dx+ I F{u)dx>\VL\ Ynin F{s) > -00. (3.8) 

n J\u\>N2 J\u\<N2 \s\<N2 

E{u) can be written in the form: 

Eiu) = ^\\ufy + J'{u) (3.9) 

with 

J^{u) = / F{u)dx. (3.10) 

It follows that E{u) is bounded from below on V, namely, 

Eiu) > ^\\u\\l + Cf, (3.11) 

where Cj := \Q\ min|s|<Af2 F{s). It's easy to see that constant Cf depends only on / and 
f2. Therefore, there is a minimizing sequence Un & V such that 

E{un) -> miE{u). (3.12) 

It follows from (13.111) that m„ is bounded in V. It turns out from the weak compactness 
that there is a subsequence, still denoted by m„, such that m„ weakly converges to v in 
V. Thus, V ^ V. We infer from the Sobolev imbedding theorem that the imbedding 

V C H^{Q) ^ L^{Q) ( 1 < 7 < ^^1) is compact. As a result, m„ strongly converges to 

V in L'^{fl). It turns out from the assumption (F2) that T{un) -^(^^)- Since \\u\\y is 
weakly lower semi-continuous, it follows from (13.121) that E{v) = mfE{u). 

The proof is completed. □ 
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4 Extended Lojasiewicz-Simon Inequality 

In what follows, we prove a suitable version of extended Lojasiewicz-Simon inequality 
required in the proof of our main result. 

Let ip he a critical point of E{u). We consider the following linearized operator 

L{v)h = -Ah + f{v + ilj)h (4.1) 

with the domain being defined as follows. 

Dom{L{v)) = {/i e H'^{n) n H^{T) : -Ay/i + d^h + h |r= 0} V. (4.2) 

The equivalent norm on V is 

\\u\\v := ||ii||i?2(n) + ||ii||if2(r). (4.3) 

It's obvious that V C L'^{fl) is dense in L'^{fl), and L{v) maps V into L^(0). In analogy 
to [29, Lemma 2.3], we know that L{v) is self-adjoint. 

Associated with 1/(0), we define the bilinear form b{wi,W2) on V as follows. 

b{wi,W2) = / (Vwi • Vw2 + f {ip)wiW2)dx + / (V||Wi • V||W2 + W1W2) dS (4.4) 

Then, the same as for the usual second order elliptic operator, L(0) + XI with A > being 
sufficiently large is invertible and its inverse is compact in L'^{Q). It turns out from the 
Fredholm alternative theorem that Ker{L{0)) is finite-dimensional. It is well known that 

Ran{L{0)) = (Xer(L(0))*)^. (4.5) 

Thus, we infer from the fact that L(0) is a self-adjoint operator that 

Ran{L{Q)) = {Ker{L{0)))^, Ran{L{Q)) © Ker{L{Q)) = L\n). (4.6) 

Next we introduce two orthogonal projections and Uji in L'^{fl), namely, n^^- is the 
projection onto the kernel of L(0) while Ur is the projection onto the range of L{0). Then 
we have the following result. 

Lemma 4.1. For 

= fR 

with Jr e L'^i^), there exists a unique solution wr G V and the following estimate holds: 

\\wr\\v <C\\fR\\. (4.7) 
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Proof. By the Fredholm alternative theorem and the regularity theorem for the elliptic 
operator (see [17]), we have a function w eV such that L(0)w = fR. Moreover w is unique 
if we require w G {KerL{Q))-^ , and (14 .7^ follows from the elliptic regularity theory. □ 

Let C{y) : V ^ L^(f2) be defined as follows: 

C{v)w = Ukw + L{v)w. (4.8) 

Then it follows from the above lemma that £(0) is bijective and its inverse C^^{0) is a 
bounded linear operator from L^(fi) to T>. 

Lemma 4.2. There exists a small constant P < 1 depending on ip such that for any 
V eV, \\v\\H^(^n) < 13 and f G L'^{^), 

C{v)w = f (4.9) 

admits a unique solution w such that w G V and the following estimate holds, 

\\wh<C\\f\\. (4.10) 

Proof. It follows from the above lemma that £(0) is bijective and its inverse C^^{0) is a 
bounded linear operator from L^(f2) to V. We rewrite (14. 9p into the following form: 

iC-\omv) - cm + I)w = C-\0)f. (4.11) 

From the definition, we have {C{v) — C{0))w = {f'{v + 11)) — f'{ip))w. 

We infer from Sobolev imbedding theorem that for any ||f ||i^2 < [3 ^ 1, there holds 

+ <C^ll^^lk2(f^)lkb- (4.12) 

Therefore, it follows that when [3 is sufficiently small, £~^(0)(£(t') — £(0)) is a contraction 
from V to V: 

\\C-'mC{v)-C{mWv,v)<\. (4.13) 

By the contraction mapping theorem, (14. lip is uniquely solvable which implies that when 
||'i^||//2(Q) < /5, C{y) is invertible, and (14.101) holds. 

Thus, the lemma is proved. □ 

Let ip he critical point of E{u). Denote u = v + ip and 



8{v) = E{u) = E{v + ip). 
11 



(4.14) 



Let 

M{v) = -A{v + i)) + f{v + i)). (4.15) 

Then for any v eV, M{v) E L^{n). 

First, we prove the following Lojasiewicz-Simon inequality for the homogeneous bound- 
ary condition corresponding to the nonhomogeneous one (II. 4p . 

Lemma 4.3. Let ip be a critical point of E{u). There exist constants 9* G (0, |) and 
P* G (0,/?) depending on ip such that for any w eV, if < [3* , there holds 

\\M{w)\\ > \S{w) - E{^)\^~'* . (4.16) 

Proof. Let J\f : V L^(fi) be the nonlinear operator defined as follows 

^^{w) = Ukw + M{w). (4.17) 

Then J\f{w) is differentiable and 

DU{w)h = C{w)h. (4.18) 

By the result in [18], we know that 

Lemma 4.4. The mapping L°°{Q) 3 u ^ f{u) G L°°(fi) is analytic. 

It easily follows from Lemma 12.11 Sobolev imbedding theorem and above lemma that 
N{w) is analytic. Since £(0) is invertible, by the abstract implicit function theorem (for 
the analytic version see e.g. [34, Corallary 4.37, p. 172]), there exist neighborhoods of the 
origin Vri(O) C P, VF2(0) C L^(fi) and an analytic inverse mapping ^ oi N such that 
^ : Vr2(0) ^ Wi{Q) is 1-1 and onto. Besides, 

M{:^{g))=g V(7GW^2(0), (4.19) 
^(7V'(f)) =v Vf G Vri(O), (4.20) 

and in analogy to the argument in [27, Lemma 1, pp.75] (see also [10, Lemma 5.4]) we 
can show that 

\\m{g^)-m{g2)\\v <C||(7i-^72|| V(7i, (72 e 1^2(0), (4.21) 
\\U{vi)-U{v2)\\ <C\\v^-V2\\v Vi;i,t;2 G W^i(O). (4.22) 

Let 01, ...,0m be the orthogonal unit vectors spanning Ker{L{Q)). 
Since \l/ is analytic, it turns out that 

m ■■= £ {f^^i^<p}j^ (4.23) 
12 



is analytic with respect to ^ = (^i,..-,^m) with |^| sufficiently small such that Ukw = 

m 

Y,^^<P^ e 1^2(0). 

i=l 

With the aid of r(^) which is an analytic function defined in M™, we are able to apply 
the Lojasiewicz inequality. By the standard argument (see e.g. [10]), we can show that, 
there exist constants 9* G (0, |) and P* G (0,/5) depending on such that for any w eV 
with \\w\\ti < P*, there holds 

llMHIl > \£{w)-m\'-'\ (4.24) 

which is exactly (14.161) . The details are omitted. □ 

Now we are in a position to prove the following extended Lojasiewicz-Simon inequality 
with boundary term. 

Lemma 4.5. Let ip be a critical point of E{u). Then there exist constants 6 G (0, |), 
/3o G (0,/5) depending on ip such that for any u G H^iVt), if \\u — tPWh'^q.) < Po, the 
following inequality holds, 

||M(i;)|| + II - A||M + 9,u + u|U2(r) > \E{u) - E{iIj)\^-' . (4.25) 

Proof. For any u G H^{fl), let v = u — ip. Then v G H^{Q). 
We consider the following elliptic boundary value problem: 



-Aw = —Av, X E Q, 

-A\\w + d^w + w = 0, X G r. 



(4.26) 



Since Av G L^{^), similar to the previous discussion for L{0), it follows that equation 
(I4.26P admits a unique solution w E V. From the if^-regularity for (I4.26P (see e.g., [17, 
Appendix Lemma A.l]), it turns out that 

||w||H2(f7) + lkllH2(r) < C\\Av\\ < C\\v\\H2(n). (4.27) 

Hence, there exists /? G (0,/3) such that for ||f ||//2(n) < /? we have 

\\w\\v<P*- (4.28) 

Here /3* is the constant in Lemma [43l Thus, (14.16p holds for w. 
On the other hand, (14.261) can be rewritten in the following form 



-A{w -v)=0, X eVI, 

- A|| {w — v) + di,{w — v) + {w — v) = A\\v — dyV — v, a; G F. 
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(4.29) 



Again from [17, Appendix Lemma A.l], there holds 



\\w - vWhhq) + \\w - v\\Hi{r) < C\\Aiiv - d^v - v\\l2^y)- 

By straightforward computation, 

\\M{w)\\ < {\\M{v)\\+C\\v-w\\h^^;) 

< {\\Miv)\\+C\\A^\v-d,v-v\\L^r)). 

Meanwhile, it follows from Newton-Leibniz formula that 



(4.30) 



(4.31) 



< 



E{w + ip) - E{v + i)) 



M{y + t{w — v)){v — w)dxdt 




'0 Jn 
•1 




(1 - t){A\\v - duv -v){v- w)dSdt 







< C {\\M{v)\\ + \\A\\v - d,v - v\\mr)) \\A\\v - d,v - t;|U2(r) 



Since 



> \E{v + ij)- E{iIj) l^-^" - I E{w + ^^J)- E{v + ^j) \^-^ 
and 0<e* <l, 2(1 - 0*) - 1 > 0, then we infer from ^M>-^M> that 



C{\\M{v)\\ + ||A||t; - d,v - v\\mr)) > \E{u) - E{^)\ 
Taking e e (0, 9*) and Po G (0,/5), such that for ||f 11^2 < Po, 

^\E{v + ^)-Em''>^- 
Let 9 = 9* - e e (0, |), then for ||f ||//2 < Pq, there holds 

\\M{v)\\ + ||Ap - d,v - v\\L2^r) > \E{u) - E{^)\'- 

which is exactly (14.251) by the definition of v. 



(4.32) 



(4.33) 



(4.34) 

(4.35) 
□ 



5 Convergence to equilibrium and convergence rate 

After the previous preparations, we now proceed to finish the proof of Theorem ll.li 
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Part I. Convergence to Equilibrium 

From the previous results, there exists an increasing sequence {tnjneN, tn — > +00 and 
e uj{uo) such that 

hm \\u(tn,x) - ij{x)\\H3in) = 0. (5.1) 

On the other hand, it follows from (12.90 that E{u) is decreasing in time. We now consider 
all possibilities. 

(1) . If there is a to > such that at this time E{u) = E^ip), then for all t > tg, we 
deduce from (12.91) that = 0. On the other hand, it follows from (II. ip (11.31) 
that 

-Ajj, = ~ut, X eCl, 

dy^ + yU = -Ut, X eV. 

Then by (12.71 ) we have 

Iktllx = / l^utdx+ / fiUtdS = \\iJ,\\ffi^^y (5.3) 
Jn Jt 

This implies that \\ut\\x = 0, i.e., u is independent oft for all t > to. Since u{x,tn) ip, 
then ( fTTBll holds. 

(2) . If for all t > 0, E{u) > E{ip), and there is to > such that for all t > to, 
V = u — Ip satisfies the condition of Lemma HTSl i.e., ||n — iPWh^^q) < Po, then for the 
constant 9 G (0, |) in Lemma [475l we have 

- |(E(..) - EWf = -e{E{u) - EWf-'^. (5.4) 



From (O), M{v) = /i. Then it follows from ( fTTiD (CT and Lemma [45] that 

- -(Eiu) - E{i,)f > g ",, II H ^ CeMHH^y (5.5) 

11/^11 + llA*llL2(r) 

Integrating from to to t, 



(£;(«) - E{ij)Y + Ce [ M\HHn)dr < (^(^(to)) - ^(^))^ (5.6) 
Since, E{u{t)) - E(^/') > 0, we have 

/ \\fJ'\\m{n)dT < +00, V t > to- 

Jtn 



(5.7) 



'to 

Thus, (Q dSITI) imply that for all t > to 



ll-Utllx^T < +CXD, (5. 
to 
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which easily yields that as t ^ +00, u(t,x) converges in X. Since the orbit is compact in 
H^{Q) n if^(r), we can deduce from uniqueness of limit that Ol.lGp holds. 

(3). It follows from (15. ip that for any e G (0,/?o), there exists N & N such that when 

n> N, 

\\ u{tn, ■) - i' \\x < \\u{tn,-) - ij\\H3(n) < ^, (5.9) 
^{E{u{t^)) - < (5.10) 

Define 

t„ = sup{ t>tn \ II m(s, ■) - ^ \h^{Q)< Po, Vs G [t„,t]}. (5.11) 

(15. ip and continuity of the orbit in H'^{Q) yield that t„ > t„ for all n > N. 
Then there are two possibilities: 

(i) . If there exists uq > N such that = +00, then from the previous discussions in (1) 
and (2), ffllBll holds. 

(ii) Otherwise, for all n > A^, we have tn < in < +00, and for all t G [tn,^n], E{ip) < 
E{u(t)). Then from (15.61) with to being replaced by and t being replaced by in we 
deduce that _ 

r \\utUdr< Cg{E{u{tn)) - E{i,)Y < |. (5.12) 

Thus we have 

||M(^n) - ^ llx < \\u{tn) - ip\\x + \\ut\\x dr < E, (5.13) 
which implies that when n +00, 

uiin) Ip in X. 

Since IJt>5^('^) relatively compact in H'^{Q), there exists a subsequence of {u(in)}, still 
denoted by {u(in)} converging to in H'^{fl)- Namely, when n is sufficiently large, we 
have 

||M(^n) - ^||//2(n) < Po, 
which contradicts the definition of t„ that \\u(in,-) — i'\\H^(n) = Po- 

Part II. Convergence Rate 

For t > toi it follows from Lemma HTSl and (15.51) that 

^{E{u) - E{^)) + C {E{u) - E{^)f^-'^ < 0. (5.14) 
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As a result, 

E{u{t))-E{^)<C{l + t)-^/^^-^'\ Vt>to. (5.15) 
Integrate (|5.5lj on {t, oo), where t >tQ, then we have 

WutWxdT <C{l + t)-'^^^-^'\ (5.16) 

By adjusting the constant C properly, we obtain 

Mt) - ^llx < C(l + t)-'/^'-''\ t > 0. (5.17) 

Based on this convergence rate we are able to get the same estimate for convergence rate 
in higher order norm by energy estimates and proper differential inequalities. 
Next we proceed to estimate ||m — -i/^Hy. 

It follows from (ll.ip - (11.4l) and the stationary problem (11.181) that 

^(u — ih) = Au, 

dV ^' ^' (5.18) 

/^ = -A(«-V^) + /(«)-/(V'), 



with the boundary condition 



-A||(m -'?/')+ dy{u - ip) + {u - ip) = n 
{u - il))t + dy^ + /i = 0. 



(5.19) 



Using (l5.18p ( l5TT9l ). we take the inner product in 7i of A ^{u — ip)t with {u — ip) to obtain 

Ij^W^- + II V« - V^f + / (/(«) - fimu - ^)dx 



|2 



+ ||V||(m - V')|L2(r) + ll'« - ^llL2(r) 
= 0. (5.20) 

On the other hand, by (15.18l) ( l5TT9ll and taking the inner product in 7i of (m — ip)t with /x, 
we have 



^QllVw-Wf + j F{u)dx- j f{%l))udx+]^\\V\\{u-il))\ 



2 

L2(r) 



1„ ^,,2 



+ 2!!^ - ^llL2(r) j + l|V/i|| + ||/i||L2(r) 
= 0. (5.21) 

Adding dOOl ) ( lOTll together, we have 

^ Qlk - ^111 + ^11^ - ^WIht) + 1 ||V||(m - + ^11 Viz - vv^f 
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+ / F{u)dx- / F{ip)dx+ / f{i))i)dx- / f{ip)udx 



+ ||V^i-V^f + ||V||(7i-^)||^,(r) 
{f{u)-fm{n-i^)dx. 



\u 



i2(r) + ||V/if + 



(r) 



(5.22) 



In what follows, we shall use the uniform bounds obtained in Lemma ETH Without loss 
of generality, we set 5 = 1 in Lemma r2.1[ 
The Newton-Leibniz formula 

F{u) = F{il)) + f{i)){u-i))+ f [ f{szu + {l~sz)ij){u-ijYdsdz, (5.23) 

Jo Jo 



yields that 



F(u)dx — / F('ip)dx + / f{il))il)dx — / f{'ip)udx 




1 /•! 



n Jo JO 



f'{szu + (1 — sz)^){u — ^Ydsdzdx 



< max \\f'(szu + (1 — sz 

s,2G[0,l] 

< C(||Vn- V^||||n-^|| -f 

< h\Vu-V^f + C\\u-' 



\U 



t > 1. 



(5.24) 



and in a similar way, we have 



(f(u)-f{mu-^)dx 



f'{su + (1 — s)ip){u — ipydsdx 




n Jo 



< -\\Vu-Vijf + C\\u 



t > 1. 



(5.25) 



Let 



\u 



11 + ^11^ - ^Whr) + \ ||V||(m - ^)||',(r^ + ^||Vn - V^f 



(I5.24p indicates that there exist constants Ci,C2 > such that 



+ / F{u)dx- / F{i))dx+ / f{ip)i)dx- / f{ip)udx 



(5.26) 



yi{t)>Ci\\u-^fy-C2\\u 



t > 1. 



(5.27) 



On the other hand, 



\u — 



< C\\u-ij\\v\\u- 
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< e\\u-i,\\l + C{e)\\u-ij\\l. (5.28) 

From f lCT7ll(l02l )f l5:M ) fl5:Wll(l5:2^ after taking e > sufficiently small, we can see that 
there exists a constant 7 > such that 

j^yi{t) + 72/i(t) <C\\u- i^Wl < C(l + t > 1. (5.29) 

As a result, 

yiit) < yi(l)e^(i-*)+Ce-^*^\l + r)-2^/(i-2e)^^ 

Jo 

< Ce~''' + Ce~^' |^^%^^(l + r)-2^/(i-2^)dr + y^*e^-(l + r)-2^/(i-2^)rfr 

< Ce--" + Ce'^' (^ei' (1 + t)-^'/^'-^'Ut + 0(1 + t)''m-2e)^jt 

< C{l + t)-^'/^^-^'\ t>l. (5.30) 
fOTll (15:281) f CTl) imply that 

< y,(t) + C2e\\u~^\\l + C2Cie)\\u-^\\i. (5.31) 
Taking e > sufficiently small, it follows from (15.17^(15.301) that 

\\u-^\\v <C{l + t)-'/^^~^'\ t>l. (5.32) 

By the C°° regularity of the solution, we are able to get the estimate for convergence rate 
in higher order norm. 

Differentiating (ll.ip -( fL4l) respect to time t respectively, we have 

utt = Afit, X efl, (5.33) 

= -^ut + f\u)ut, xeQ, (5.34) 

Utt + d„Ht + Ht = 0, X G r, (5.35) 

fj^t = —^\\Ut + di,Ut + Ut, X G r. (5.36) 
Multiplying ( I5.34p by Ut and integrating by parts on Q, using ( I1.3p ( l5.36l l. we get 

2jt ("^'""' + w^'WW)) + 11^^*11' + (I vii^*r + ^0 ds 

19 



f'{u)uidx. (5.37) 
Assumption (F3) yields that there is a certain positive constant Mf > 1 such that 

f'i'S) > -Mf, s e M. (5.38) 

Thus, 

- [ f{u)uldx < M/||utf . (5.39) 

It follows from (II. ip that 



— / Vut -V ^dx — I ^UtdS - ""• 
< \\^Ut\\\\Vii\\ + \\fi\\mr)\\ut\\mr) 



L2(r) 



< e\\Vutf + e\\ut\\h^^^ + ^11 V/if + ^ll/^ll W), (5.40) 

In (EHOD, taking 
it follows that 

~ (llV/if + ll^lli.(r)) + l\\Vu,r + II V||t.,||i.(r) + i|k,||i.(r) 

< C^(l|V/xf + ||/i||i.(r))- (5.42) 

Multiplying (15.341) by Uu and integrating by parts on Q, using (15.33^(15.351) (15. 36p . we get 

+ ll/i*llW) + l|V/i.||^ 

f"{u)v!ldx. (5.43) 



By (15.401) and Lemma [2Tn the righthand side of (15.431) can be estimated as follows 

f"{u)u\dx 



n 



< Ci\u\L^)\\ut\\l, < C\\u\\h2 {\\Vutp\\u,p + \\Ut 

< l\\vu,r+c\\u,r+c\\u,r 

< l\\vu,r + c\\u,r 



< 

- 4 



\ + |k,||i.(r)) + C|| V/zf + Cll^lli.(r), t > 1. (5.44) 
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Then (15.430 becomes 



1 d 
2dt 



||Vn,f + / /'(«Krfa;+||V||n,||i.(r) + ||n,||i.(r) 
Jn 

< ^(l|Vn,f + h,||i.(r))+C||V/if + C||/.||i.(p), t>l. (5.45) 
Multiplying fl5.45l) by Ei G (0, 1] and adding the resultant to ()5.42p . we obtain 



2 

L2(r) 



+^ilkJi2(r)) + ^ (llVwif + ||Mt||i.(r)) + W^wm^^D 
+^i||/^dli2(r) + ^i||V/iiir 
< C*(||V/.f + t>l. (5.46) 



Let 



?/2(t) = ||V/i||'+||^||i2(r)+5i||VMt|p+£i^/'(M)M2rfx+£i ||V||Mt||^,(r)+£i||Mi||i.(r) (5.47) 
It follows from Lemma 12.11 that 

y2{t) <C, t>l. (5.48) 

Taking 

^1 = ^, (5.49) 

we can deduce from ( 15.39p -f l5.4ip that 

y2{t) > I (llV/if + Ml.^r)) + ei (^||Vn,|r + || V||«,||^,(,) + ^h*lli.(r)) • (5.50) 
Now we take k > such that 

<1 + C*)<1. (5.51) 

Next, we multiply (15:461 ) by n and add the resultant to (15:221 ). then ( iOOll (lOHD ( 15:321 ) 
yield that there exists a constant 7 > such that 

j^Mt) + Ky,{t)] + 7[l/i(t) + ^y2m <C\\u- < C(l + t)-2^/(i-2^), t > 1. (5.52) 

Similar to (15.301) . we have 

2/1 (t) + ^^y2{t) < C(l + t)-^'/('-^<^\ t > 1. (5.53) 
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Hence, from (15:271) fOOl) (021) (l533l) we know 

y2{t)<C{l+t)-^"^^-^'\ t>l, (5.54) 
which together with (I5.50p gives the following 

M\mm + \\ut\\m(^n) + \\ut\\m(^T)<C{l + t)-"^^-^'\ t>l. (5.55) 
By the elliptic estimate (see [17, Corollary A.l]), 

\\u - '^Wn^n) + 11^ - '^Whht) < C + ||/(^^) - fi^)\\m{n) + WfJ'WmiT)) ■ (5.56) 

Lemma [2.1! and Sobolev imbedding theorem imply that 

\\f{u) - fmimm < C\\^i - ^llv, t > 1. (5.57) 

On the other hand, from (15 ■21) . the elliptic regularity theory and Sobolev imbedding 
theorem, we have 

WfJ'WnHT) < C\\fi\\H2(n) < C (^\\ut\\ + \\ut\\^i^^^^ < C\\ut\\min)- (5.58) 
As a result, we can conclude from (I5.32p . (I5.55I) - (I5.58I) that 

\\u - nn^n) + 11^ - nH^T) < C(l + t)-'/^'-''\ t > 1. (5.59) 
Summing up, the proof of theorem 11.1! is completed. 

Remark 5.1. Following the same method, we can continue to get estimates of convergence 
rate in higher order norm. 

Remark 5.2. We notice that, in order to get the convergence rate estimates ( I5.32p ( I5.55p 
(I5.59p . we have to use the uniform bound for the solution in higher order norm, e.g. 
Lemma \2.1\ which is not valid for t = 0. Thus, the constant C in (I5.32p(j5. 550 (15.590 
depends on 6 in Lemma \2.1\ More precisely, for any 5 > w;e have 

\\u - + \\u - + WutWv < Cs{l + t)-"^^-^'\ yt>S. (5.60) 

the constant Cs depends on \\uo\\v and 6. Moreover, 

lim Cs = +00. (5.61) 
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